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Abstract. We present an exactly solvable toy model which describes the emergence of a pseudogap in an 
electronic system due to a fluctuating off-diagonal order parameter. In one dimension our model reduces 
to the fluctuating gap model (FGM) with a gap A(x) that is constrained to be of the form A{x) = Ae l ® x , 
where A and Q are random variables. The FGM was introduced by Lee, Rice and Anderson [Phys. Rev. Lett. 
31, 462 (1973)] to study fluctuation effects in Peierls chains. We show that their perturbative results for 
the average density of states are exact for our toy model if we assume a Lorentzian probability distribution 
for Q and ignore amplitude fluctuations. More generally, choosing the probability distributions of A and Q 
such that the average of A(x) vanishes and its covariance is (A(x)A*(x')) = A 2 S exp[— \x — x'\/£\, we study 
the combined effect of phase and amplitude fluctutations on the low-energy properties of Peierls chains. 
We explicitly calculate the average density of states, the localization length, the average single-particle 
Green's function, and the real part of the average conductivity. In our model phase fluctuations generate 
delocalized states at the Fermi energy, which give rise to a finite Drude peak in the conductivity. We also 
find that the interplay between phase and amplitude fluctuations leads to a weak logarithmic singulatity 
in the single-particle spectral function at the bare quasi-particle energies. In higher dimensions our model 
might be relevant to describe the pseudogap state in the underdoped cuprate superconductors. 

PACS. 71.23.-k Electronic structure of disordered solids - 02.50.Ey Stochastic processes - 71.10.Pm 
Fermions in reduced dimensions 

1 Introduction controversial. Several mechanisms have been proposed. Ac- 

cording to Schmalian et al. [Q the normal state of the 
The physical origin of the pseudogap behavior observed in underdoped cuprates can be modeled by a nearly antifer- 
the normal state of the high-temperature cuprates is still romagntic Fermi liquid, and the experimentally observed 
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pseudogap behavior is closely related to strong antifer- = Z\^e _ ' r_r . (4) 



romagnetic spin fluctuations. An alternative explanation, 
which has been advanced by Emery and Kivelson re- 
lates the pseudogap behavior to precursor superconduct- 
ing fluctuations. In this scenario thermal fluctuations of 
the phase of the superconducting order parameter are re- 
sponsible for a destruction of superconductivity above the 
transition temperature T c . However, in a wide range of 
temperatures T > T c the local amplitude of the super- 
conducting gap is finite. In this paper we shall propose 
a simple exactly solvable phcnomcnological model which 
describes the destruction of phase coherence due to phase 
and amplitude fluctuations of the superconducting order 
parameter in the pseudogap state. 

To study superconducting fluctuations in a normal metal 
one can start with the Gorkov equation for the 2x2 ma- 
trix Green's function for electrons with energy dispersion 
e(k) that are coupled to a space-dependent complex pair- 
ing field A{t) §, 

[u - H r ] g( d =V ( r , r',u) = 8(r - r> , (1) 

/ \ [w - H x ] Q(x, x',u) = S(x - x')a , (5) 

e(-iV r )-fi A(r) \ 
H r = (2) 

\ A*(r) e(iV r )-M/ . -iv F d x A(x) \ 

H *=\ \ ■ ( 6 ) 

Here, Co is the 2x2 unit matrix and /i is the chemical I ^\* ( x \ i Vf Q x J 

potential. In the absence of true superconducting long- We shall refer to Eq ^ ^ as the H amiltonian of the fluctu- 
range order the pairing field A(r) can be considered as a ating gap model ( FGM ). A11 quantities depend now para- 
random variable with zero average and correlations that metrically on r± and n physical observables should be 
fall off exponentially with distance, averaged over all directions of n. In this paper we shall 
(Z\(r)) = (3) oru y consider the effective one-dimensional problem de- 
fined by Eqs. (ra) and (ph. We require that the first and 

(A(r)A* <r')) ^ f n^ 3 ^ A(r)A* (r>) 

\ K ' v J T>{A}e~ s { A } the second moments of the fluctuating gap A(x) are given 



Here, is the Ginzburg-Landau functional of the or- 

der parameter field, £ is the correlation length, and the 
energy scale A s characterizes the strength of the correla- 
tions. 

To simplify the algebra and to make contact with other 
theoretical work on pseudogap physics, we shall focus in 
this work on the semiclassical limit of the Gorkov equa- 
tion, which are related to the so-called Andreev equation 
[Q. In the weak coupling limit, where |^i(r)| is small com- 
pared with the chemical potential, we may linearize the 
energy dispersion in Eq. (Q) for wave-vectors k close to 
the Fermi surface, provided we are only interested in long- 
wavelength, low-energy properties of the system. In the 
semiclassical limit it is useful to decompose the position 
vector as r = xn + r± where n is a unit vector in the 
direction of the momentum of the electron, and rj_ is or- 
thogonal to n. Writing d x — n • V r , Eqs. (|l|) and (^) can 
be replaced by an effective one-dimensional problem || 
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by One easily verifies that Eqs. ([?]) and (||) are indeed satis- 



fied. Note that for commensurate Peierls chains the order 
parameter field can be chosen real, while it is complex 
for incommensurate chains. In this work we shall focus on 
In the following, we shall construct a special non-Gaussian the incommensurate case, where zero-energy states and 



(A(x)) = , (7) 
(A(x)A*(x')) =Ale-\*- x 'M . (8) 



probability distribution of A{x) satisfying Eqs. (0) and (||) the associated Dyson singularities are absent Lee, 

for which Eq. (||) can be solved exactly. Moreover, as will Rice and Anderson treated the effect of the order parame- 

be briefly discussed in Sec. |, it is straightforward to gen- ter fluctuations on the average electronic density of states 

eralize our model to dimensions d > 1 and to arbitrary (DOS) (p(u))) within the Born approximation. Within this 

energy dispersions e(k), although the calculation of phys- approximation one finds that, in the regime where the di- 

ical quantities becomes more tedious. mensionless parameter 

Apart from its relevance in the semiclassical theory vf 
of superconductivity, the problem defined by Eqs. (ra) to 

x-i ... . is small compared with unity, the DOS develops a pseu- 

(|8|) describes also the low-energy physics m quasi-one- 

i. nra T dogap for \cu\ ^ A S7 with a minimum given by IBI 

dimensional Peierls and spm-Pcicrls systems fg,|6||. Lee, tJ 

Rice and Anderson |j| used this model to study fluctu- (p(0)) p _ 7 ^ 

Po y/l + 7 2 ' 

ation effects close to the Peierls transition. In this case 

Here, 

A(x) can be identified with the fluctuating Peierls order , 

Po = (13) 

parameter, and the two diagonal elements in our Hamilto- F 

r,, i • ,i is the DOS for A(x) = 0, which is a constant due to the 

represent the kinetic energy of the electrons m the v ' 

.... r , n „ . . , , , . , linearization of the energy dispersion. Note that Eq. (112) 

vicinity of the two Fermi points ±kp. Physical quantities 

, , , . , , , , .... ,. , ., ,. predicts for 7 < 1 to leading order 

should again be averaged over the probability distribution 

of A(x), which can be obtained from the Ginzburg-Landau (p(0)) p ^ - ^ ^-1 

Po 

expansion ||. Within the Gaussian approximation, the 
truncated Ginzburg-Landau functional in the disordered 
phase is of the form 



which disagrees with a non-perturbative result by Sadovskii 
[[To), who found for the model defined by Eqs. (||) to (||) 
for a Gaussian distribution of A(x) 

Po 



where 



However, the algorithm constructed by Sadovskii 1 1 



A q - J dxe tqx A(x) . (10) nQt exact QHj^ so that it is not clear whether Eq. (|l|) 
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is correct or not. To clarify this point, we have recently Finally, it should be mentioned that a generalization 

developed an exact numerical algorithm for calculating of the model defined in Eqs. (||) to (|sj> has been used in 

the DOS of the FGM || . For a Gaussian distribution of Ref . Q to explain the pseudogap behavior in the cuprates 

A(x) with zero average and covariance given by Eq. (||) within antiferromagntic Fermi liquid theory. Then the scalar 

the result is field A{x) should be replaced by a matrix field J2i Si(x)&i, 

, , „ „ where <7, are the Pauli matrices, and the fields Sdx) repre- 

«a[27] b ocr b , (16) 

Po sent the components of the antiferromagnetic spin density 

n field. In fact, the recent interest in the non-perturbative 

where 

approach invented many years ago by Sadovskii [jfO| is 
a = 0.6397 ± 0.0066 , b = 0.6397 ± 0.0024 . (17) motivated by its possible relevance to the cuprate super- 
conductors. 

Hence, for Gaussian disorder with a finite correlation length 
both perturbation theory and Sadovskii's algorithm do not 

give the correct ^-dependence of the average DOS at the 2 ExaCt Green's function of the fluctuating 

Fermi energy. Another attempt to investigate the discrep- gap model for A(x) = Ae l ® x 
ancy between Eqs. ( p2| ) and ( |T5| ) numerically was recently 

In this section we shall solve Eq. (|J) exactly for a special 

made by Millis and Monien |Q. They found for the ex- 

form of the probability distribution of A(x) which is con- 

ponent b in Eq. ( |16|) a value between 2/3 and 1, which 

structcd such that its covariance is given by Eq. (||). To 

is outside our error-bars in Eq. (p"7j). Note, however, that 

begin with, let us perform the following gauge transfor- 

Millis and Monien studied a lattice regularization of the 

mation [13], 

continuum model (|6|), and no attempt was made to care- 
fully relate the bare parameters that appear in the lattice Q(x, x' , uj) = e^ a ^ a3 Q(x, x' , ui)e~% a ( x , (18) 
and the continuum models. In this work we shall show 

where the gauge function a(x) will be specified shortly. 

that the exponent characterizing the behavior of the DOS 

From Eq. (§) we find that the transformed Green's func- 

at the Fermi energy on £ is non-universal in the sense that 

tion Q(x,x',ui) satisfies 

it depends on the precise form of the probability distribu- 

r vp da{x) -ia(x) 
tion of the fluctuating gap. In particular, the non-Gaussian u> h ivpOxVz — A{x)e y 'a + 

terms in the Ginzburg-Landau functional can change the _ A* (x)e ia ^ Q(x, x' , u>) = S(x — x')<7q . (19) 

numerical value of this exponent, so that the behavior 

Suppose now that A(x) is of the form 

given in Eqs. ([l6]) and ( |l7| ) can only be expected to be 

correct for Gaussian disorder. A(x) = Ae 1 ^ , (20) 
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where A and Q are both random but independent of x. (^) and (Q) if A and Q are interpreted as random vari- 

Then the ^-dependence of A[x) in Eq. ( p^) can be re- ables. To obtain the exponential decay of the covariance 

moved by choosing a(x) = Qx. Moreover, with this choice we require that the probability distribution of the random 

the second term on the left-hand side of Eq. (J19|) reduces momentum Q is a Lorentzian, 



to a constant 



ro — — 



1 



vp da(x) v F Q 
~2 dx~ _ 2 



(27) 



(21) 



so that 



[u> — r\ + ivpd x a^ — Aa + — A*<j-] Q(x, x' , to) 

= 8(x-x')a ■ (22) w i t h 

Thus, a phase of the order-parameter varying linearly in 



or equivalently for the random energy shift rj defined in 
Eq. ©, 

„ 7 1 



7r ry 2 + 7 2 



2£ 



(28) 



(29) 



space can be absorbed by a finite shift of the energy. Eq. The random variab i e A shomd be distributed such that 
(E2) is translational invariant and is easily solved by a 



Fourier transformation, 



(A) a = , 
{\A\ 2 )a = A 2 s , 



(30) 
(31) 




(23) 



(24) 



where (■ ■ -)a denotes averaging over the probability dis- 
tribution of A. From Eqs. ( p7[ ) to ( pl| ) it is then easy to 
show that the first two moments of the distribution of 
A(x) are indeed given by Eqs. (|?|) and (|8|). Note that Eqs. 



Combining Eqs. (yj), Q23|) and (|2J) and defining 

g(q,q',uj) = JdxJ dx'e- l ^ x -«' x ')g(x,x',u) . 
we finally obtain 



(25) 



G(q,q',u) = 



/ 2i:8(q - q')[m - 2?7 + v F q] _ 

lu — 2rj + VFq] [oj — v F q] — \A\ 2 [lu — 2rj + vpq] [lu — vpq] — \A\ 2 



0J) and (31) include the cases of pure phase and pure am- 
plitude fluctuations. To describe pure phase fluctuations 
we choose A = A s e llp , where the phase if is uniformly 
distributed in the interval [0, 2tt). Then 

27r<5(<7 -q' - Q)A 



\ 



2tt<5(<7 - q' + Q)A* 2ir6{q -q')[uj-2r)- v F q] 

V [lu — 2rj — v F q][uu + v F q] — \ A\ 2 [lu — 2r\ — v F q][uu + v F q] — \ A\ 2 I 

(26) 



The crucial observation is now that, in spite of the sim- 



ple form (20) of A(x), it is still possible to satisfy Eqs. 



I A 



dip 
2^ 



(32) 
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Since physical quantities should be independent of the bution of 77 is a Lorentzian of width 7 is therefore not in 

constant phase ip and therefore should only depend on contradiction to the work of Rcf . ||l6f . Obviously, our pa- 

\A\, the process of averaging amounts to replacing \A\ by rameter 7 in Eq. ( p9| ) is the analog of the parameter W in- 

A 8 . To take into account amplitude fluctuations we follow troduced in Eq. (9) of Ref. Note, however, that Franz 

Sadovskii [fi~4[[l(| and choose a Gaussian distribution for and Millis [[l6) did not consider amplitude fluctuations of 

the real and imaginary parts of A, the order parameter, which are described by our second 



\ am 



dReA dlmA -\a\ 2 /a 2 random variable A. As noted above, Gaussian amplitude 

1 I *A\ e ■ (33) 

fluctuations with a probability distribution given by Eq. 

The disorder averaging of any functional T{A(x)\ is de- ,rr^. , , ... . „ . . .. 

( p3| ) have been studied many years ago by Sadovskii [I14J. 

fined by Thus, in the present work we combine the models intro- 

(F{A{x)}) = (^J dQVg F{Ae tQx }^ . (34) duced by Sadovskii @ and by Franz and Millis (T§ such 

that we take both amplitude and phase fluctuations into 

What is the physical meaning of an order parameter of the 

account and still obtain an exactly solvable model. 

form (pp|)? In a superconductor such an order parameter 

In the following section we shall calculate a number of 



describes a state with a uniform superflow 15 . The gauge 

physical quantities for this model exactly and confirm the 

transformation ( |l8| ) corresponds to choosing a coordinate 

intuitive picture J^,0 that phase fluctuations fill in the 

system where the superflow vanishes; r\ is the associated 

gap at the Fermi energy and render the system metallic. 

energy shift. A more detailed physical justification for such 
a spatially constant random energy shift r\ in the normal 

state of the cuprate superconductors has been given by 3 Calculation of physical quantities 

Franz and Millis [|l6) : they pointed out that within a semi- 

3.1 Single-particle Green's function and spectral 

classical approximation the effect of the quasi-static fluc- 

function 

tuations of the phase of the order parameter field A(x) can 

be described by such an energy shift n. Franz and Millis _ , .. r zFTTK , 1 , ,1 rr 

J 6J ' Because (A) = 0, it follows from Eq. (|26|) that the off- 

lllql also presented a perturbative calculation of the prob- ,. , , r,,-, ,r 

4 — ' diagonal elements of the disorder averaged Green s func- 

ability distribution V r , of 77, using earlier results by Emery . . , , , , . . 

' tion vanish, and that the diagonal elements are 

and Kivelson Qj. Because in Ref. [|l6| a cumulant expan- 
sion of Vri was truncated at the second order, the form of (Gaa{q, q , w)) — 2nS(q q )G a {q,io) , (35) 
V v was found to be Gaussian by construction. However, 

where 

there are certainly non-Gaussian corrections to the form 

r , \ / g - 2?7 + av F q \ 

of V v given in Ref. |lq| . Our assumption that the distri- "W'^ \ r w _ 2rj + avFq][u - avFq] - \A\ 2 / ' 
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Here, a = + refers to Q\\, and a — — refers to The 

averaging over the Lorentzian distribution ( |28| ) of the ran- w — av Fl 



dt ^ . (40) 



dom energy shift rj can be performed analytically, Recently Kuchinskii and Sadovskii @ arrived precisely at 

G a (q,LO + l0+) = 



1 \ Eq. ( |40| ) within a diagrammatic attempt to estimate the 

= \ |^4p / ' 

\ui — avpq — — — — qj^^vp A accuracy of the method developed in Ref. |l(| for Gaus- 

(37) 

sian disorder. For a better comparision with Sadovskii's 

where (. . .) A denotes averaging over the probability dis- Green , s function ca i cu i ate d in Ref. @, let us represent 

tribution of A. In the case of pure phase fluctuations, as Eq Q as a continued f ract i n. 
described by Eq. (p3), this averaging is trivial, so that 



1 

ui — avpq — Sa a (q, oj + i0 + ) 



G£(q,w + iQ + ) = — TT , (38) 



with the self-energy given by 

Sg h (g, " + = - A2 °. VF - (39) 
uj + avpq + i-g- 

Eq. (|9|) agrees precisely with the lowest order Born ap- 
proximation, which was used in the seminal work by Lee, 
Rice, and Anderson We have thus found a special 
probability distribution of A(x) where the lowest order 
Born approximation for the average single-particle Green's 
function is exact: the order parameter is in this case of the 
form A(x) = A s e %( ^ x+%<p , where Q has a Lorentzian dis- 
tribution of width l/£, and the random phase ip merely 
assures (A(x)) = 0, but due to gauge invariance does not 
affect any physical quantities. 

On the other hand, if in addition to phase fluctua- 
tions also amplitude fluctuations are important, there are 
corrections to the Born approximation. For Gaussian am- 
plitude fluctuations given by Eq. (^) we find after sub- 
stituting t = \A\ 2 /A 2 S 

G^ +am (q,UJ + lO+) = 
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Expressing the integral on the right-hand side of Eq. ( pio| ) 
in terms of the incomplete r~ function and using the known 
continued fraction expansion of this function p8[ |, we ob- 
tain for the self-energy 



oj + avFq + i 



^ 



oj — avpq 



oj + avFq + i—r 



.v F 2Aj 



oj - av F q jj^- 

oj + avpq + i— 



(41) 



For the same model with Gaussian disorder the algorithm 
due to Sadovskii p~0| ] produces the continued fraction ex- 
pansion 

Sl lldovskii (q,LJ + iO + ) = 



.Vf 

oj + avpq + i— 



oj — avFq + 2i 



£, vf 2A 
oj + avFq + 3i— 



oj — avpq + 4i— r 



5 w + avFq + 5i— ... 

(42) 

Note that only the first two lines in Eqs. ([Id]) and 
agree. Kuchinskii and Sadovskii argue in Ref. jl7| that 
the true behavior of the Green's function for Gaussian 
disorder lies somewhat in between Eqs. pi] ) and piq). In 
our model, the coexistence of amplitude fluctuations with 
phase fluctuations (which are related to our random en- 
ergy shift rj) generates a completely new feature in the 
average spectral function. The latter is related to the av- 
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erage Green's function via tron. Note that such a singularity is weaker than the alge- 

2 braic singularities that are typically found in the spectral 

2w6(q - q')(p(ak F + q,u)) = Im (G aa (q, q' ,u + i0 + )) , 

7T 

function of a Luttinger liquid. Of course, such a weak sin- 

(43) 

gularity cannot be called a quasi-particle peak. It is impor- 

Using Eq. @ we find 

tant to point out that in the presence of amplitude fluctu- 
(p(akF + q, to)) p + = ations alone or phase fluctuations alone such a logarithmic 

/ dt- ^tt^ , (44) singularity does not exist. Recall that for pure phase flue- 
Jo (t - to 1 + q z ) z + A-y z (LU - aq) z 

tuations our model has the same spectral function as pre- 
dicted by the Born approximation for the self-energy 
while for pure amplitude fluctuations our model reduces 
to the model discussed by Sadovskii in Ref. |l4|]. Note 
also that the approximate spectral function produced by 



it A s J Q (t - u> 2 + q 2 ) 2 + 47 2 (d> - acq) 

where q = vpq/A s , to — uj/A s , and 7 = Vf/(2A s £). Rep- 
resentative results for different values of 7 are shown in 
Figs, [j] and |[ The dashed line is the spectral function for 
7 = (i.e. without phase fluctuations), which is easily 
calculated analytically, 



Sadovskii's algorithm [|19|,|20| for Gaussian disorder with a 
(p(akp + q,to)) = finite correlation length does not exhibit any logarithmic 

A~ 1 0(lo 2 — q 2 )\& + ag|e - ^ 2 ~* 2 ) . (45) singularities. Whether an exact calculation of the spectral 

function for more realistic probability distributions could 

The important point is now that for any finite 7 the spec- 
confirm this result or not remains an open question. 

tral function exhibits a logarithmic singularity at u) = 

From Fig. || it is clear that the line-shape of the spec- 

avFq- In the vicinity of this singularity the leading behav- 

tral function in the vicinity of the singularity is rather 

ior of the spectral function can be calculated analytically. 

broad and asymmetric. Such a behavior has recently been 

In the regime 

seen in the photoemission spectra of a one-dimensional 



' A 2 £ A 2 

\w - av F q\ < min <{ —f-, —5 \ (46) band-insulator 01 

vf \uj + avFq\ 



the integral in Eq. (|44|) can be approximated by 

(p(rf P +^)) ph+ "~ -^-hJ 1 
nA* 



2j\lo — aq\ 



VF dn 



3.2 Average density of states 



The average DOS is defined by 



(47) 



vf\oj — av F q\ 

Thus, the interplay between phase fluctuations (described 
by our random phase factor e z< ^ x ) and amplitude fluctua- 
tions (described by random fluctuations of \A\) gives rise Tr {G{x, x, ui + i0 + )) . . 

v F \y/\A\ 2 - {iu + i-f) 2 I A 

to a logarithmic singularity at the bare energy of the elec- (49) 



(p(w)) = --ImTr{g(x,x,u } + iO+)) . (48) 



Performing the g-integration in Eq. (|37| ) we find 

1 / LU + 47 
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where 7 is given in Eq. and \fz denotes the principal A numerical evaluation of R(j) is shown in Fig. ^ For 
branch of the square root, with the cut at the negative small and large 7 we obtain to leading order 

y/n^f , 7 < 1 
1 , 7> 1 

absence of amplitude fluctuations (see Eq. (|32j)) we may For large £ the DOS at the Fermi energy is 



real axis. Note that phase fluctuations simply generate an 
imaginary shift ij to the frequency in Eq. (|l9ft. In the 



(56) 



replace \A\ — > A s in Eq. ((49|), so that we obtain for the Mj m 



average DOS 



((0(0))P ™ „ _y_ j ^ >>Zis; (57) 
which should be compared with the result obtained within 



= Im 7= , , (50) 



Po \J\ — z 2 the Born approximation, see Eq. (]14|), 

where we have defined = ^ ^ 1 

w + 17 _ ._ . 
Z A s 10 l ^ ' Hence, Gaussian amplitude fluctuations increase the value 

Eq. agrees exactly with the perturbative result by of the DOS at the Fermi energy as compared with pure 

Lee, Rice, and Anderson For u = we recover Eq. @. phase fluctuations. However, from Fig. § it is evident that 

On the other hand, in the presence of additional Gaussian the qualitative behavior of the DOS is correctly predicted 

amplitude fluctuations, with probability distribution given by a model with pure phase fluctuations, which exactly 

by Eq. (|33|), we obtain reproduces the perturbative result ||. Let us emphasize 

(n(uj)) ph+a - m f°° e~ l z that this is not the case if A(x) has a Gaussian distribu- 

= Im / * * . (52) V ' 

tion: the prediction of lowest order perturbation theory, 
(p(0)} oc is in disagreement with the exact numeri- 
cal result for Gaussian disorder, (p(0)) oc £ -0 - 64 (see Eq. 



Po Jo Vt 



z 



A numerical evaluation of Eq. (p52) is shown in Fig. a. For . 
7 = the integral in Eq. (|5^) can be done analytically ] , ( /• r 1 



and reduces to the result obtained by Sadovskii Il4), which r— K 

(|16|)). We thus conclude that the behavior of the average 

does not contain phase fluctuations. In this case the DOS 

DOS at the Fermi energy of the FGM in one dimension 

vanishes quadratically for small frequencies, . . . . 

is non-universal and sensitive to the detailed form of the 



\am 



2lo , < 1 . (53) probability distribution of A{x). 

3.3 Lyapunov exponent and localization length 



Po 

For any finite £ the DOS at the Fermi energy (i.e. at uj = 0) 
is finite. From Eq. (p2[) we find 



( n(()))P h+a - m Since the energy dispersion of the FGM is linear, the 

A-^L =i?(7), (54) 

/3 ° Schrodinger equation H x ijj w {x) — ujipufa) is a system of 

with . . 

linear first order differential equations. Fixing the two- 

r°° e-* 

R(fY) = 7 / dt - . (55) f. . . . . , . 4 .. , 

Jq _|_ 72 component wave-function ^(x) arbitrarily at one space 
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point xq therefore constitutes the wave-function at all A plot of the inverse localization length £~ 1 (u>) ph+a - m is 
points x. In a disordered system, the Lyapunov exponent given in Fig. ^[ For 7 — > only amplitude fluctuations are 
k(oj) characterizes the exponential growth of the magni- left, and Eq. (j63J) reduces to 



tude of the wave-function at large distances \x — Xq\ J22| , 
\ipu>(x)\ ~ \ip u (xo)\exp[n(cj)\x - x \] . (59) 



e , 7 — > . (64) 



A s £(uj) am 2 
In the presence of phase and amplitude fluctuations the 
Strictly speaking, the Lyapunov exponent is defined by the S eneral ex P ression © simplifies at the Fermi energy to 



limit \x — xq\ — > 00 of this equation and assumes a cer- ^£ _ p(=\ fce^i 

ZU(0)P h + am ~ K1) ' K ' 



where the dimensionless function P("y) is given by 



tain value with probability one p2| . In one dimension the 
inverse of the Lyapunov exponent can be identified with 

/>oo 

the mean localization length. According to the Thouless P(l) = / dte" 1 \/t + j 2 — 7 . (66) 



formula the mean localization length £(u) can be obtained 
from the real part of the disorder-averaged single-particle 
Green's function. Originally the Thouless formula was de- 
rived for a one-band model with quadratic energy disper- 

' \/tt/2 , 7 < 1 

sion j23| , but it can be shown to hold also for the FGM, P(j) ~ i . (67) 



A comparison of Eq. ( pq ) with the corresponding expres- 
sion obtained from Eq. ( |62| ) for phase fluctuations is shown 
in Fig. H For small and large 7 the leading behavior is 



l/(2 7 ) , 7 » 1 
In the white noise limit £ — * 0, A s — > 00 with A 2 £ = const 
~du JhJ) = ReTr ^( a; ' x ' w + * 0+ )) • ( 60 ) only the behavior of P(7) for large 7 matters, and in this 

limit both Eq. ( p2[ ) and Eq. ( |65| ) reduce to the known 
white-noise result 



where it can be written as [p4[|2 
_d 1 

Integrating the Thouless formula for Eq. (E^), we obtain 



Re(V|A| 2 -(c + z7)2)^-7, (61) 



^ " "' U " UF ^ wrth ^ = const . (68) 



where the constant of integration is uniquely determined 
by the requirement lim^^oo £~ 1 (to) = 0. For pure phase 
fluctuations Eq. (|6l]) reduces to 



^(0) 2 7 v F 
An extrapolation of this white-noise result towards finite 
correlation lengths is shown as the dotted line in Fig. 
[| Evidently, for large 7 the behavior of the localization 



I'F 



A s e(ui)p h + a 



j ^ )p . — Re \/\ — (<D + Z7) 2 — 7 , (62) length becomes independent of the precise form of the 

probability distribution of the disorder. For 7 1 the lo- 
calization length begins to deviate significantly from the 
7 . white-noise limit and approaches a finite value of the order 
(63) of vf/A s for 7^0, the precise value of which depends 



while with additional Gaussian amplitude fluctuations 

I'F 



Re 



dte~ l y/t - (u> + 27) 2 
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on the type of the disorder. We emphasize that for a real polarization (n(q,iuj m )), which is given by 

order parameter the low-frequency behavior of the local- \ f dp f dp' 

2 7 rd(q-q J ){n(q,iu> rn ))=-~J2J —J — 

ization length is dominated by the Dyson singularity, so ™ 

xTr(g(p + q,p' + q',iQ n+m )Q(p',p,iuj n )) . (69) 

that in this case 1/^(0) = for any finite value of 7, see 

Refs. [p^,|2q] . Here, (3 is the inverse temperature, uj m = 2wm/f3 are 

To compare the localization length of our exactly solv- bosonic Matsubara frequencies and u n = 2n{n + \)/(3 

able toy model with phase and amplitude fluctuations are fermionic ones. Given the average polarization, the 

with the case where the distribution of A(x) is a Gaussian, average conductivity is easyly obtained from 
we have evaluated the Thouless formula (pOf) numerically , . ,. iu> „,., 

for Gaussian colored noise with correlation length £, us- 

j— 1 . . In this work we shall only consider the real part of the 

mg an algorithm 129] similar to the one developed m Kef. 

in. . . , . ... , conductivity at q = 0, 

|pfl. I he numerical results for vf/(A s i{0)) are shown as 

the open circles in Fig. |. In view of the simplicity of our Re / a (u)) = lim Re (a(q, u)) = e 2 uj lim ( Inl ^(g'^ + ?0+ )) 

g^O q— >0 q^ 

model the agreement with Eq. (pq) is quite spectacular. 



Hence, the localization length of our model with phase and Substituting Eq. @ into Eq. © and performing the 
amplitude fluctuations is a very accurate approximation Matsubara sum, we obtain for the average polarization 
to the localization length of the FGM with Gaussian dis- 



(II(q,iu m )) = 

dp E p E p+q + £ p € P+q + \A\ 2 



order. The dashed line in Fig. [sj describes the localization 

length for the case where we ignore amplitude fluctuations \ J 2tt 2E p E p+q 

'f(E p -y)-f(E p+q - V ) 



in our model, which is equivalent to the perturbative re- 
sult by Lee, Rice, and Anderson ||. The agreement with fi^p + 7 ?) ~ fj^p+R ^ 

E p — E p+q + iLo m 

the case of Gaussian disorder is not so good, in particular p ^p E p E p+q — i p £, p + q — \A\ 2 



E p — E p+q — iuj r , 



in the pseudogap regime 7 ;$ 1. J 27r 2E p E p+q 



l-f(Ep-r))-f(E p+q + r]) 



E p + E p+q — iL0 ri 



\- f(E p + r,)- f(E p+q -r,) 



(72) 



3.4 Average conductivity E p + E p+i + luJm 

where we use the notation E p — + l^l 2 ) 1 / 2 , £ p = v F p 
The DOS and the spectral function [see Eqs. © and @] and f(E) = l/[eJ m + 1] is the Fermi-Dirac function. Set- 
involve only the diagonal elements of the single-particle ting r\ — in Eq. (|7^) we recover Eq. (2.10) of Ref. Q. 
Green's function. The simplest physical quantity which Expanding Eq. (|7^) for small q and performing the aver- 
involves also the off-diagonal elements of Q is the average age over the Lorentzian distribution of 77, we obtain in the 
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limit of zero temperature (/? — > oo), found that the localization length £(0) at zero frequency is 

ne 2 7r / , \ finite. In fact, we believe that for Gaussian disorder with 

Re(<7(o;)> = — -(V|A| 2 +7 2 -|A|) 5(u>) 

, . moments given by Eqs. (m) and ml) the conductivity of 

I ne2 axctanf M ) \ A \ 2 °^-\ £1 \ ,( 73) 

m V 7 / \ cj2 \fu 2 — \A\ 2 j A ' the one-dimensional FGM does not exhibit a Drude peak, 

where n/m = vp /n and 7 is defined in Eq. (E9J). For pure because the eigenstates at u> = should all be localized 



phase fluctuations the averaging over the distribution of for a given realization of the disorder [g2|,|l£j|. On the other 
A is trivial and simply leads to the replacement \A\ — > hand, for our choice A(x) = Ae l( ^ x with spatially constant 
A s . Then the conductivity exhibits a Drude peak with but random A and Q, the Green's function is not self- 



weight given by j~ 1 (y / A 2 + 7 2 — Z\ s ), which is separated averging, so that its spatial average is not identical with its 

from a continuum at higher frequencies by a finite gap disorder average. As a consequence, there is a finite prob- 

A s . Gaussian amplitude fluctuations wash out the gap ability of finding delocalized states at the Fermi energy: 

but do not remove the Drude peak. Averaging over the for \w — 77] > |^4| the solutions of the Schrodinger equation 

probability distribution of the amplitude A given in Eq. are simply plane waves, whereas for \u — rj\ < \A\ there 

( |33| ) we obtain is a gap in the spectrum, and the Schrodinger equation 



Re(cr(w)) 



9 

ne 



in 



7rD(7)*M + ^-C(7,w) 



does not have any normalizable solutions. Hence, depend- 

(74) 



ing on the realization of the disorder, the system is either 

where we have used again the notation 7 = 7 / A s , uj = . , . , , „ . ^ 

a perfect conductor or an insulator. .Because m Eq. (|6 

lo/A s , and the dimensionless functions D(-j) and C(7,(D) 



we have defined the inverse localization length in terms of 

1 f°° 

D(j) = - dte-*[y/t + i* - Vi] , (75) 
7 Jo 



the disorder averaged Green's function, the value of £ 1 (uj) 



is determined by those realizations of the disorder where 



C(7 ui) = arctan [ ) \uj\ [ dte"^ 1 — - — (76) localized states at energy uj do not exist. However, the 
V 7 / Jo VI ~t 

A graph of D(fy) is shown in Fig. [?]. Physically D(ff) is the 



dimensionless renormalization factor for the weight of the 
Drude peak, with D = 1 corresponding to an unrenormal- 
ized Drude peak. The leading terms in the expansion of 
D(7) for small and large 7 are 

^r7 , 7 « 1 

' { 2 ■ (77) 

1 , 7> 1 



probability of finding delocalized states at the Fermi en- 
ergy is finite, and can be expressed in terms of the function 
P(7) defined in Eq. (p 



Wdeloc(O) = (0( V 2 - \A\ 2 )) 

= 1 - -?=P(7) 

\/7T 



At the first sight the existence of a Drude peak in our J ^7 i 7^1 



(78) 



model is rather surprising because in Sec. 3.3 we have 1 , 7 > 1 



14 Lorenz Bartosch and Peter Kopietz: Exactly solvable toy model for the pseudogap state 

A graph if Wdcioc(O) is shown as the dashed-dotted line order the real part of the conductivity is known to vanish 
in Fig. [?]. Note that the qualitative behavior of Wdcioc(O) for small frequencies as uj 2 ln 2 (l/u;) || . Thus, apart from 
is very similar to the weight D of the Drude peak. The the logarithmic correction, the incoherent part of the con- 
conductivity of quasi-one-dimensional Peierls systems be- ductivity of our simple model shows the generic behavior 
low the Peierls transition (for which (A(x)) 7^ 0) has been of one-dimensional disordered electrons. Note also that for 
discussed in Refs. [^6[^7). The authors pointed out that small 7 the relative weight of the Drude peak is of the or- 
in this case a gapless collective mode associated with flue- der of 7, so that the incoherent contribution dominates, 
tuations of the phase of the order parameter generates a The w hite-noise limit is defined by letting A„£ -> 
finite Drude peak. In our toy model, r, describes such a while keeping A*£ finite. In this case D(fy) approaches 
gapless mode. unity. In fact, in the white-noise limit the average conduc- 

As discussed in Sec. g our model is also relevant to tivity is not modificd by the d i sor der at all because the 
describe higher-dimensional systems such as superconduc- f unct i n Z\j 1 C( 7 , Q) vanishes if we let A s -> 00. 
tors within a quasiclassical approximation. In this case it 
is physically reasonable to expect that phase fluctuations 
of the superconducting order parameter generate delocal- 4 Conclusions 
ized states at the Fermi energy P,[l6|. Then we indeed 
expect a finite Drude peak in the conductivity, which is 

In this work we have introduced a simple exactly solvable 

broadened by disorder and becomes a sharp ^-function in 

toy model which describes the combined effects of phase 

the superconducting state. 

and amplitude fluctuations of an off-diagonal order pa- 
Let us now focus on the incoherent part of the conduc- 

rameter on the physical properties of an electronic system. 

tivity, which is described by the dimensionless function 

Although we have only discussed the one-dimensional ver- 

C(7, uj) in Eq. (|74|). A graph of this function is shown in 

sion of this model with linearized energy dispersion, the 

Fig. ||. For large correlation lengths, i.e. 7 <§; 1 there are 

exact solubility of our model does not depend on these 

three characteristic regimes where C("f, uj) can be approx- 

features, so that our calculations can be generalized to 

imated by 

more realistic models of electrons in dimensions d > 1 



C(7,0) 



3 7 uj , \uj\ <C 7 with non-linear energy dispersions. In this case the fluctu- 



47T 

6 



'-\uj\ , | 7 | < \u\ < 1 ■ (79) ating gap should bg chosen f the from ^( r ) _ Ae iQ r . To 

f \uj\~ 3 , 1< |w| satisfy ( A ( r y = o and (A{r)A*{r')) = A 2 s e-\ r - r '\^ , the 

For 7 <C \uj\ this agrees with the result of Ref. [^4). Note random variable A should be distributed such that Eqs. 

that for a one-band model with Gaussian white noise dis- (^FJ) and ( ^lf ) are satisfied, while the distribution Pq of 



the other hand, even very weak three-dimensional interac- 
tions can lead to a phase transition leading to long-range 
order and a finite Drude peak as found in our model. We 
expect that forward scattering by disorder (which we have 
ignored in our calculation) will broaden the Drude peak 
[ [30| . Experimentally, peak structures in the far infrared 
well below the pseudogap regime have been observed in 
the optical conductivity of several quasi one-dimensional 
Peierls systems above the Peierls transition [Bl] . 
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the d-dimensional random-vector Q should be strictly one-dimensional disordered system, the disorder 

,„ _ 1 / , -iQ-r -|r|/£ / sft \ should lead to the localization of all eigenstates, resulting 
Q _ (2n) d J 6 ' 

in a vanishing zero temperature dc conductivity jp9|. On 

For d = 1 this reduces to Eq. @, but in d > 1 Eq. @ 
is not a Lorentzian. 

In one dimension our model describes the disordered 
phase of Peierls and spin-Pcicrls chains. We have presented 
explicit results for the density of states, the localization 
length, the single-particle spectral function, and the real 
part of the conductivity. Let us emphasize three points: 

(a) The mean localization length of our toy model, 
which we have defined via the Thouless formula (|6^), is an 
excellent approximation to the mean localization length of 
the FGM with Gaussian disorder. Although the respective Our model also describes superconducting fluctuations 
density of states agree quite well on a qualitative level, de- in d > 1 within a scmiclassical approximation. Recall 
viations become substantial for large correlation lengths, that our Eq. (J^) for the Green's function in d = 1 is 
leading to a different scaling behavior as a function of £. formally equivalent to the Andreev equation for the semi- 

(b) The interplay between phase and amplitude flue- classical wave-function of a superconductor. The latter can 
tuations gives rise to a weak logarithmic singularity in the be obtained from the more general Gorkov equation (|l]) 
single-particle spectral function of our model. Whether in the limit of a slowly varying order parameter. To cal- 
this singularity is just an artifact of our toy model or not culate physical observables, the solutions of the Andreev 
remains an open question. equations should be averaged over the classical trajecto- 

(c) The conductivity of our model exhibits not only ries of the electrons Q], which we have not done in this 
a pseudogap below the energy scale A s but also a Drude work. Therefore we cannot make any quantitative corn- 
peak at oj = with a weight that vanishes as 1 /£ for parisons with experimental data for high-temperature su- 
£ — > oo. While the qualitative picture of the continuous perconductors. However, some qualitative features of our 
part should be generic for more realistic one-dimensional results seem to agree with experiments. In particular, in 
disordered systems (up to logarithmic corrections for small our model the pseudogap in the conductivity coexists with 
frequencies p9[|), the Drude peak in our model is due to a small Drude peak. Such a behavior has been seen exper- 
the existence of delocalized states at the Fermi energy imentally in the normal state of high-temperature super- 
which are created by phase fluctuations. However, in a conductors [|32f| . In our model the Drude peak is a direct 
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Drude peak. In this respect our model describes a bad (1999). 
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Fig. 3. Average DOS <p(a;)) ph+am [see Eq. @] as a function Fig. 5. Inverse localization length r 1 (u) ph+aln /v' 1 A s [see 
of G) = u/A s for 7 = (dashed line), and j = 0.1, 0.25, 0.5, 1. Eq. (^)] as a function of uj — uj/A s for 7 = (dashed line), 



For smaller 7 the pseudogap becomes deeper. 
1 




and 7 = 0.1, 0.25, 0.5,1 
1 




Fig. 4. Solid line: numerical evaluation of R{"f) 



Fig. 6. Inverse localization length P(j) = Vf/A s £(Q) at the 



(p(0)) ph+am /Po as a function of 7 = vf / (2A S £) (see Eq. (p5[)). Fermi energy for different types of disorder. Solid line: phase 

Dashed-dotted line: the same quantity without amplitude flue- and amplitude fluctuations, see Eqs. ( |65|) and (|66|); dashed- 

tuations (see Eq. (^)), which amounts to calculating the av- dotted line: phase fluctuations, see Eq. (p^); dashed line: ex- 

erage DOS from the self-energy in first order Born approxima- trapolation of white noise limit P(j) = 1/(27), see Eq. (^)- 

tion, as was done by Lee, Rice, and Anderson 0|. The circles The circles are numerical results for Gaussian disorder, ob- 

are numerical results for Gaussian disorder, obtained via the tained via an exact numerical algorithm J25[|. 



exact numerical algorithm of Ref. 
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Fig. 7. Solid line: dimensionless renormalization factor D of 
the weight of the Drude peak as a function of 7 = vf/(2A s £), 
see Eq. (j7jj); dashed-dotted line: probability Wddoc(O) for find- 
ing delocalized states at the Fermi energy, see Eq. (TTsb - 




12 3 4 



Fig. 8. Incoherent part C(7,oj) of the conductivity as a func- 
tion of G) = ui/A a , see Eq. (|t^) . From top to bottom: 7 = 
(dashed line) and 7 = 0.1, 0.25, 0.5, 1. 



